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SECTION  1 
INTRODUCTION 


In  order  to  characterize  a  distant  seismic  pulse  according  to  the  energy 
which  generated  it,  a  source  function^  is  needed  which  can  be  used  to  initiate  the 
seismic  signal  at  sufficient  range  from  the  event  that  the  subsequent  propagation  can 
be  described  in  linear  terms  according  to  the  properties  of  the  earth.  Thus  an  ob¬ 
served  seismic  signal  implies  a  source  function  which  may  be  compared  with  known 
functions  for  discrimination  and  yield  estimation.  Therefore,  the  properties  of  such 
source  functions  must  be  known  to  use  this  technique.  Near-field  data  from  various 
events  have  been  taken.  The  smallest  strains  observed,  due  to  practical  free  field 
instrument  placement  selection,  are  typically  no  less  than  10“®  and  are  often  much 
greater.  If  the  pulse  at  this  extreme  range  undergoes  no  further  significant  nonlinear 
modification,  its  characterization  can  supply  the  needed  source  function.  However, 
if  any  additional  nonlinear  changes  are  important,  a  useful  source  function  cannot 
be  determined.  Consequently,  it  is  important  to  characterize  any  possible  nonlin¬ 
ear  attenuation  of  moderate  strain  pulses,  preferably  through  methods  which  allow 
generalization  to  all  placement  media  of  interest. 

There  is  a  long  history  of  development  of  numerical  methods  for  calculation 
of  strongly  nonlinear  behavior  in  the  near  field  of  UGTs.**’  The  gross  nonlinear  ef¬ 
fects  of  vaporization,  crushing,  cracking  and  plastic  behavior  induced  by  strains  much 
greater  than  10~^  must  be  taken  into  account  to  obtain  an  understanding  of  near  field 
data  taken  at  NTS.  Complex  equations  of  state  or  constitutive  relations  including 
effective  stress  and  porosity  are  needed  for  a  description.^  These  methods  generally 
apply  for  strains  greater  than  10~^;  at  lesser  strains  linear  behavior  is  assumed.  How¬ 
ever,  cis  indicated  in  the  following  section,  it  seems  clear  that  some  residual  nonlinear 
effects  continue  out  to  strains  less  than  10“®.  Ideally  one  would  like  to  have  a  phys¬ 
ically  based  model  of  any  nonlinear  effects  in  order  to  use  existing  near  field  data  to 
define  a  more  distant  linear  source  function  or  at  least  to  determine  what  effect,  if  any, 
the  nonlinearities  over  the  moderate  strain  range  have  on  the  linear  source  function. 

In  order  to  determine  this,  we  are  attempting  to  develop  nonlinear  consti¬ 
tutive  relations  for  the  mildly  nonlinear  attenuation  in  the  moderate  strain  regime 
(scaled  ranges  of  10*  m/kT^/*  to  10^  m/kT**^*  or  strains  from  10"*  to  10"®)  for  explo¬ 
sively  generated  seismic  pulses.  The  relations  are  to  be  used  to  determine  the  character 
of  seismic  pulses  in  the  linear  low-strain  regime  (beyond  10*  m/kT^/*)  which  may  then 
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serve  as  source  functions  for  regional  or  teleseismic  propagation.  A  linear  source  func¬ 
tion  can  be  found  by  numerical  propagation  of  an  experimentally  well  known  initial 
pulse  through  the  moderate  strain  range  to  account  for  mild  nonlinear  attenuation. 

Experimental  data  provide  significant  restrictions  on  candidate  constitutive 
relations.  Attenuation  in  a  salt  medium  frmn  Salmon,®"*®  Cowboy,**"*®  Livermore 
small  scale  explosions,*®  Rockwell  damped  oscillation  experiments*®  and  New  Eng¬ 
land  Research  ultrasonic  pulse  attenuation,*®  with  one  exception,  are  fairly  consistent 
internally.  The  first  three  provide  a  detailed  indication  of  the  change  of  amplitude 
and  sha^  of  explosively  driven  pulses.  Proposed  constitutive  relations  must  repro¬ 
duce  these  data.  The  data  from  the  wide  range  of  yields  in  salt  indicate  that  yield*^® 
scaling  applies  with  a  remarkable  precision;  if  all  times  and  distances  are  scaled  by 
yield*^®,  the  scaled  amplitudes  and  shapes  of  pulses  from  all  experiments  are  nearly 
the  same.*'  Thus  the  initial  pulse  on  entering  the  moderate  strain  regime  must  scale  at 
some  small  scaled  range  and  the  subsequent  attenuation  must  result  from  constitutive 
relations  which  have  no  time  or  space  scales  which  are  fixed  by  the  medium.  This 
provides  a  significant  reduction  in  the  domain  of  allowable  nonlinear  behavior.  For 
example,  the  relations  may  be  a  function  of  the  strain  but  the  may  not  depend  on  the 
strain  rate. 

In  order  to  investigate  nonlinear  constitutive  relations,  a  standard  numerical 
time  stepping  method  is  used.  The  technique  which  we  have  used  is  that  of  taking 
the  observed  Salmon  initial  velocity  pulse  at  small  range  (166  meters)  as  a  source 
and  comparing  the  resulting  pulses  as  they  are  propagated  through  material  subject 
to  candidate  constitutive  relations.  The  results  are  compared  with  observed  signals 
at  larger  ranges.  For  any  constitutive  relation  the  effective  Q  associated  with  the 
attenuation  may  be  determined  but  it  must  be  emphasized  that  nonlinear  attenuation 
cannot  be  properly  described  by  a  Q  function;  still  Q  may  sometimes  be  useful  for 
comparison  with  past  work.  The  fundamental  comparison  of  the  data  with  calculations 
is  not  in  terms  of  the  Q  but  in  terms  of  reproduction  of  waveform  including  both 
amplitude  and  shape. 
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SECTION  2 


MODERATE  STRAIN  ATTENUATION  DATA  IN  SALT 


The  most  comprehensive  attenuation  data,  over  the  range  of  moderate 
strains,  exists  in  the  medium  of  salt.  For  explosive  sources  the  laboratory  results 
of  Larson'^  covers  strains  from  10“^  to  10"^,  the  Salmon  field  test®  covers  strains  for 
10~®  to  10“^  while  the  Cowboy  field  test^^  series  covers  strains  from  10“^  to  10“®. 
The  Rockwell*®  laboratory  decaying  oscillations  complete  the  salt  data  by  ranging 
from  10“®  to  10“®.  The  NER*®  laboratory  ultrasonic  pulse  propagation  experiments 
go  from  less  than  10“®  to  more  than  10“®.  The  essential  results  of  these  experiments 
will  now  be  reviewed. 


Larson’s  laboratory  data  were  taken  from  the  effects  of  a  series  of  small 
chemical  explosives  in  blocks  of  pressed  salt.  Velocity  data  were  taken  at  several 
gauges  such  that  the  totality  of  the  examples  from  all  shuts  covered  peak  velocity 
to  compressional  velocities,  which  is  approximately  the  peak  strain,  of  10“*  to  10“®. 
Using  a  triplet  of  records  from  a  single  shot,  it  was  estimated  that  over  peak  strains 
of  1.4  X  10'®,  7.0  X  10~^  and  4.6  x  10“^,  the  Q  changes  from  12.5  to  24.9.  The  corner 
frequency  for  this  small  scale  experiment  was  about  5  x  10^  Hertz. 


The  Salmon  data  were  generated  by  a  5.3  kT  nuclear  explosion  in  salt.^  The 
work  by  McCartor  and  Wortman^  as  well  as  McLaughlin  and  Gupta®  shows  that  the 
high  quality  data  clearly  indicate  a  high  level  of  attenuation.  It  is  estimated  that  for 
peak  strains  from  4  x  10~®  to  3  x  10“^,  the  effective  Q  is  to  order  of  10  at  a  corner 
frequency  of  about  6  Hertz.  This  Q  appears  nearly  constant,^  perhaps  increasing 
mildly  with  range,®  over  the  order  of  magnitude  strain  range  available.  However,  in 
view  of  the  fact  that  small  strain  data  from  other  experiments,  including  the  decoupled 
event  Sterling®  which  used  the  Salmon  cavity,  indicate  a  much  larger  Q,  it  seems  likely 
that  there  are  residual  nonlinearities  at  the  extreme  range  of  the  Salmon  data  so  that 
Q  must  increase  at  larger  ranges  (and  so  smaller  strains) . 


Tittmann*®  has  studied  the  attenuation  of  flexural  and  torsional  harmonic 
oscillations  of  salt  samples.  It  has  been  found  that  the  attenuation,  expressed  in 
terms  of  Q~*  as  a  function  of  strain  amplitude,  tends  to  be  nearly  constant  for  strains 
from  10~®  to  10“®,  then  increase  for  greater  strains.  The  attenuation  is  a  decreasing 
function  of  confining  pressure  with  an  average  Q  of  approximately  200  at  10“®  strain 
at  a  frequency  of  400  Hertz  for  confinement  consistent  with  the  explosions. 
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The  regime  covered  by  the  COWBOY  experiments  overlaps  that  for  Salmon, 
by  beginning  at  strains  of  about  5  x  10“*  and  extending  out  to  10“®.  The  COWBOY 
experiments  consisted  of  a  series  of  shots  over  a  range  of  yields,  both  coupled  and 
uncoupled,  such  that  the  individual  events  generally  had  only  a  few  instruments  in 
operation.  Several  authors  have  combined  these  data  by  invoking  the  experimentally 
compelling  evidence  for  simple  scaling,  which  suggests  that  the  velocity  pulse  data 
from  various  yields  can  be  equated  by  scaling  all  times  and  distances  by  the  cube  root 
of  the  yield.  This  is  particularly  dramatic  when  peak  velocity  is  plotted  against  scaled 
range  giving  a  consistent  curve  over  about  ten  orders  of  magnitude  in  yield.  Larson 
has  pointed  out  that  the  full  velocity  pulses,  as  well  as  their  peak  values,  tend  to  be 
preserved  through  scaling. 

Trulio**’  has  analyzed  the  scaling-combined  Cowboy  data  to  estimate  Q  as 
dictated  by  decrease  in  peak  displacement  with  range  by  using  one  decade  at  a  time  in 
the  frequency  domain.  He  has  found  that  attenuation  decreases  as  range  increases  in 
a  fashion  inconsistent  with  linearity  expressed  as  dispersive  harmonic  potential  waves. 
If  the  coupled  Cowboy  data  are  fit  assuming  a  Q  independent  of  range  (a  possibility 
due  only  to  scatter  in  the  data),  Q  must  be  a  function  of  frequency  ranging  from 
about  5  at  low  frequencies  (1  Hz  at  Salmon  scale)  to  nearly  100  at  high  frequencies 
(32  Hz  at  Salmon  scale).  However,  simple  scaling  is  inconsistent  with  this  fit  since 
scaling  requires  linearity  coupled  with  a  Q  independent  of  frequency.  If  the  decay 
of  the  reduced  displacement  potential  is  fit  by  the  form  exp(-u;r/2cQ),  it  is  seen 
that  simple  scaling  can  result  so  long  as,  for  constant  c,  Q  is  a  function  of  a  variable 
which  is  unchanged  by  scaling  -  such  as  u)r.  Trulio  has  shown  that  the  a  range  of 
data  from  Salmon,  through  Cowboy  and  Cowboy  Trails,  do  have  the  property  that  an 
effective  Q  can  be  expressed  as  a  function  of  wr  for  all  the  experiments.  Obviously  this 
indicates  that  at  fixed  frequency,  the  attenuation  expressed  in  terms  of  Q,  is  dependent 
upon  range.  Since  it  is  assumed  that  the  medium  is  approximately  homogeneous, 
the  implication  is  clear  that  the  attenuation  must  be  amplitude  dependent  and  so 
nonlinear. 


Minster  and  Day^*  have  used  the  scaled  peak  velocity  data  for  COWBOY 
to  determine  if  these  data  require  and  amplitude  (nonlinear)  or  frequency  dependent 
Q  for  consistency.  They  determined  that  a  Q“*  which  consists  of  a  small  constant 
(consistent  with  small  strain  data  from  Tittmann)  plus  a  term  proportional  to  the 
peak  strain  provides  a  good  fit  to  the  data.  The  observed  attenuation  effects  do  not 
firmly  indicate  the  need  for  a  frequency  dependent  Q  but  indicate  that  an  amplitude 
dependent  Q  provides  a  more  convincing  fit  than  a  constant  Q.  Based  on  the  small 
strain  Q,  which  is  several  hundred  as  seen  in  the  Rockwell  experiments,  it  is  concluded 
that  there  must  be  nonlinear  attenuation  in  the  COWBOY  strain  regime. 
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Coyner*®  of  New  England  Research  has  carried  out  a  series  of  laboratory 
experiments  for  which  compressional  and  shear  ultrasonic  pulses  consisting  of  about 
two  cycles  at  100-200  kHz  were  propagated  through  a  sample.  Attenuations  were 
calculated  using  a  spectral  ratio  technique.  Variation  of  the  attenuation  with  peak 
strain  amplitude  and  confining  pressure  were  determined.  Experiments  were  carried 
out  with  several  materials  including  Sierra  White  Granite,  Berea  Sandstone  and  Dome 
Salt.  For  the  dome  salt  it  was  found  that  over  a  strain  range  of  5  x  10~^  to  3  x  lO-® 
and  for  a  confining  load  range  of  0.1  to  1  MPa,  the  P-wave  attenuation  is  nearly 
constant  and  can  be  described  by  a  Q  of  about  20.  There  is  no  particular  evidence  of 
nonlinearity  in  these  data  alone. 

In  summary,  the  salt  data  appear  to  support  the  hypothesis  that  explosively 
generated  pulses  encounter  nonlinear  attenuation  for  strains  much  greater  than  10“®. 
.No  detailed  knowledge  of  the  attenuation  mechanism  currently  exists  but  there  does 
appear  to  be  a  consistency  in  than  the  explosively  generated  pulses  closely  obey  simple 
scaling  which  suggests  that  the  mechanism  must  be  rate  independent.  The  effective 
Q  for  this  process  increases  with  increasing  frequency  and  increases  with  decreasing 
strain.  In  the  next  section  we  shall  attempt  to  exploit  a  feature  of  the  Salmon  near¬ 
field  pulses  which  suggests  a  physical  mechanism  which  can  account  for  the  attenuation 
data. 
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SECTION  3 


MECHANISM  FOR  SALMON  ELASTIC  PRECURSOR 


The  Salmon  data  have  a  feature  which  may  be  useful  in  understanding  some 
nonlinear  effects  of  attenuation.  Each  of  the  pulses  experimentally  observed  at  six 
sensors  at  ranges  from  166  meters  to  660  meters  exhibit  a  discontinuity  in  the  slope 
upon  the  initial  steep  rise  in  pulse  velocity,  as  seen  in  Figure  1.  This  appears  as  a 
toe-like  behavior  in  the  leading  edge  of  the  velocity  profile  which  has  been  described 
by  Ferret®  as  an  “elastic  precursor”  to  the  main  pulse.  The  absolute  amplitude  of 
the  toe  remains  approximately  constant  with  range,  at  a  particle  velocity  of  about 
0.5  m/s,  while  its  amplitude  relative  to  the  peak  increases  with  range.  This  precursor 
amplitude  corresponds  to  a  compressional  strain  level  of  t  «  10“^.  The  leading  edge 
of  the  pulses  (i.e.,  that  disturbance  earliest  in  time)  propagates  at  a  speed  of  about 

4.7  km/sec  while  the  pulse  peaks,  always  after  the  toe,  propagate  at  a  speed  of  about 

3.7  km/sec.  The  elastic  compressional  speed  of  mild  disturbances  in  this  salt  medium 
found  from  independent  measurements  is  typically  about  4.6  km/sec.  This  indicates 
that  the  precursor  signal  seen  in  the  Salmon  data  is  due  to  elastic  behavior  while  the 
subsequent  pulse  suffers  a  lower  propagation  speed  due  to  some  relaxation  or  plastic 
behavior.  If  propagation  were  purely  plastic  the  sequence  of  pulses  seen  at  the  Salmon 
sensor  sites  would  be  as  shown  in  Figure  2. 

Ferret  suggests  an  elastic-plastic  material  behavior  might  account  for  the 
data  in  perhaps  one  of  two  ways.  First,  the  precursor  could  develop  at  large  strain, 
where  an  elastic  limit  is  exceeded  from  radii  much  smaller  than  instrumented  for 
Salmon,  and  continue  to  propagate  in  front  of  a  following  plastic  wave.  Second,  it 
may  be  that  the  precursor  develops  in  the  moderate  strain  region  if  dome  salt  has  an 
elasto-plastic  nature  at  such  strains.  In  either  case,  the  modulus  of  salt  must  be  a 
function  of  the  strain  -  that  is,  the  medium  is  nonlinear. 

If  the  precursor  develops  at  large  strain  there  must  be  an  elastic  limit  beyond 
which  plastic  behavior  provides  a  lower  modulus.  When  such  a  medium  is  dynamically 
loaded  beyond  the  elastic  limit,  a  leading  pulse  at  the  elastic  limiting  stress  is  generated 
followed  by  a  larger  amplitude  but  slower  plastic  wave.  If  the  elastic-plastic  transition 
is  not  sharply  defined,  the  resulting  pulse  could  consist  of  a  gently  rising  leading 
elastic  front  which  smoothly  merges  with  the  main  plastic  pulse  much  as  seen  in 
the  Salmon  data.  Ferret  points  out  that  if  some  energy  from  the  plastic  component 
is  fed  to  the  elastic  portion  during  propagation,  the  amplitude  of  the  elastic  piece 
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Figure  1.  Salmon  radial  velocity  pulses  seen  at  166,  225,  276,  318,  402 
and  660  meters  range  (time  zero  on  plot  corresponds  to  initial 
signal  arrival  at  166  meters  which  was  at  0.036  seconds  after 
detonation). 
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Calculated  pulses  at  Salmon  ranges  for  pure  plastic  medium 


could  remain  nearly  constant,  but  this  is  quite  speculative.  What  is  known  is  that 
a  variety  of  laboratory  experiments  show  that  strong  impulsive  sources  can  produce 
elastic  precursors  in  a  variety  of  media.  For  example,  work  of  Ahrens  and  Duvall** 
with  planar  pulses  in  quartz  generated  by  explosives  exhibits  an  apparent  elastic  stress 
limit  of  about  70  kbar  corresponding  to  strains  of  about  10“ *.  This  produces  a  leading 
edge,  described  as  the  elastic  shock,  which  propagates  at  a  speed  in  excess  of  that  of 
the  deformational  portion  of  the  pulse  which  follows.  The  elastic  shock  propagates 
with  an  equivalent  modulus  which  is  greater  that  the  static  modulus  at  this  high  stress. 
It  is  speculated  that  the  elastic  wave  is  supported  by  a  higher  than  equilibrium  shear 
stress.  After  the  elastic  component  has  passed,  the  shear  stress  is  apparently  reduced 
to  the  static  value  by  a  plastic  or  fracture  process.  This  experiment  as  well  as  others*'' 
which  show  an  elastic  precursor  consistent  provide  elastic  limits  of  many  or  tens  of 
kilobars  in  contrast  with  the  Salmon  data  which  gives  a  precursor  amplitude  of  about 
5  bars.  Consequently  this  mechanism  does  not  seem  a  likely  means  of  accounting  for 
the  Salmon  data  which  give  a  nearly  constant  and  small  precursor  amplitude  which 
seems  to  begin  near  the  166  meter  sensor  range  rather  than  be  well  developed  by  this 
time. 


The  second  possibility  indicated  by  Ferret  is  that  of  having  the  precursor 
develop  locally  in  the  observation  region  based  on  moderate  strain  plastic  behavior. 
W'hile  there  is  no  accepted  dynamical  equation  of  state  for  dome  salt,  Ferret  points 
out  that  dome  salt  is  known  to  be  highly  plastic  -  under  static  conditions,  it  is  nearly 
hydrostatic.  Thus  plastic  deformation  of  salt  at  moderate  stresses  is  apparently  nor¬ 
mal.  In  order  to  account  for  the  precursor  data  in  Salmon,  the  equation  of  state  would 
have  to  provide  linear  behavior  up  to  a  threshold  (a  threshold  of  about  5  bars,  much 
less  that  the  70  kbars  found  for  the  elastic  shock  discussed  above)  and  through  some 
deformation  of  the  material,  relax  the  modulus  abruptly  (on  the  Salmon  time  scale) 
to  a  value  which  provides  a  compressional  propagation  speed  about  20%  less  than  for 
infinitesimal  strains  in  undisturbed  material.  Having  a  threshold  at  the  observed  and 
nearly  constant  precursor  amplitude  avoids  having  to  account  for  the  constancy  of 
the  precursor  by  arguments  of  convenience  about  feeding  energy  back  from  the  plastic 
portion  of  the  wave. 

As  an  equation  of  state  model  which  can  be  consistent  with  the  Salmon  data, 
consider  a  medium  for  which  the  shear  modulus  permanently  (that  is,  does  not  recover 
until  the  pulse  is  past)  decreases  upon  having  a  critical  strain  threshold  exceeded;  the 
compressional  modulus  before  and  after  exceeding  the  strain  threshold  reflects  the 
compressional  speeds  at  the  beginning  and  peak  of  the  Salmon  pulses,  respectively. 
This  will  be  referred  to  as  a  shear  failure  model.  Depending  upon  the  relation  between 
the  compressional  and  shear  moduli,  complete  shear  failure  may  occur,  meaning  that 
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the  elastic  shear  modulus,  //,  goes  to  zero.  For  the  example  used  in  this  discussion,  the 
compressional  speed  decreases  to  80%  of  its  original  value.  We  have  taken  the  Lame 
constants  A  and  p  to  have  a  ratio  of  2.  Thus  a  decrease  of  the  compressional  speed, 
((A  +  where  p  is  the  density,  of  20%  corresponds  to  reducing  n  to  about 

38%  of  its  elastic  value  when  A  is  held  fixed.  Note  that  the  reduction  of  modulus  at  a 
fixed  strain  is  consistent  with  scaling  since  the  strain  is  a  unitless  quantity  and  there 
is  no  rate  dependence.  The  scaling  restriction  does  require  that  the  relaxation  time 
of  the  modulus  change  be  short  compared  with  any  representative  time  scale  of  the 
data;  we  take  the  transition  to  be  instantaneous.  Figure  3  indicates  the  stress-strain 
curve  which  results  from  this  model  for  the  modulus.  In  order  to  determine  the  effect 
on  the  pulse  propagation  we  used  the  observed  Salmon  velocity  at  166  meters  as  the 
source.  The  calculations  were  carried  out  using  a  standard  finite  difference  methods 
whose  details  are  provided  in  Appendix  A. 

As  an  initial  effort,  the  elastic  threshold  was  taken  at  a  compressional  strain 
of  10~‘*;  the  resulting  pulse  sequence  at  the  ranges  to  observation  stations  for  Salmon 
is  as  shown  in  Figure  4.  Note  that  the  character  of  the  calculated  precursor  is  much 
like  that  seen  experimentally,  in  Figure  1,  in  that  the  leading  feature  is  drawn  out,  the 
transition  to  the  main  pulse  takes  place  at  a  constant  amplitude  and  the  peak  now 
moves  at  a  significantly  lower  speed.  Still  the  amplitude  of  the  main  peak  does  not 
decrease  as  quickly  as  the  data  indicate. 

When  the  modulus  decreases,  the  elastic  energy  in  the  pulse  also  decreases 
in  a  manner  approximately  proportional  to  the  square  of  the  compressional  wave 
speed.  Since  the  modulus  reduction  is  permanent,  this  energy  is  lost  to  the  pulse  and 
goes  into  heating  the  medium.  For  the  parameter  used,  over  a  full  cycle  for  which 
most  of  the  pulse  exceeds  the  critical  strain,  approximately  one-third  of  the  original 
elastic  energy  will  be  lost.  This  corresponds  to  an  effective  Q  of  about  13  for  peak 
strains  well  in  excess  of  10“^  (for  small  strains  less  than  this  threshold,  there  will 
be  no  loss).  This  value  of  Q  is  far  less  than  that  expected  for  very  small  strains 
but  it  is  still  more  than  the  5  to  10  seen  for  Salmon  attenuation.  The  addition  of  a 
moderate  level  of  linear  attenuation  consistent  with  that  seen  for  small  to  moderate 
strains  in  other  experiments  will  improve  the  agreement  (for  example,  the  NER  data 
suggest  Q  ss  20).  More  importantly,  the  use  of  a  partial  shear  failure  mechanism  will 
automatically  terminate  once  the  pulse  weakens  so  that  the  peak  strain  falls  below 
the  critical  strain  threshold  value.  This  will  produce  a  sharply  changing  effective  Q  in 
a  manner  suggested  by  the  Cowboy  data  as  given  in  Figure  5.  Thus  this  first  attempt 
to  reproduce  the  Salmon  data  is  encouraging  but  there  remain  differences  which  may 
be  accounted  for  by  refining  the  details  of  the  attenuation  mechanism. 
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Figure  3.  Stress-strain  curve  for  partial  shear  failure. 
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SALMON  FROM  QF  166-660 
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Calculated  pulses  at  Salmon  ranges  for  shear  failure. 
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Figure  5.  Cowboy  Q  estimates  as  a  function  of  peak  strain.  The 
horizontal  bars  represent  the  two  record  strain  values  while 
the  vertical  bars  indicate  the  spread  of  Q  over  a  decade  of 
frequency.  The  fit  of  Minster  and  Day  is  provided.  (From  Ref. 
3). 


The  attenuation  from  partial  shear  failure  does  not  produce  an  amplitude 
for  the  pulse  at  660  meters  which  is  as  small  as  that  seen  experimentally.  More  atten¬ 
uation  can  be  added  to  attempt  to  match  better  the  data  by  employing  a  linear  Q  of 
sufficient  value.  A  method  of  inclusion  of  a  linear  absorption  band  Q  in  time  stepping 
finite  difference  methods  has  been  demonstrated  by  use  of  Fade  approximants.  Our 
application  of  this  method  is  outlined  in  Appendix  B.  This  formalism  was  employed 
using  a  target  Q  of  10  with  a  range  of  half  amplitude  frequencies  of  1  to  100  Hertz  (Q 
rises  above  20  beyond  these  values).  The  sequence  of  pulses  which  result  using  both 
the  partial  shear  failure  and  a  linear  Q  of  10,  starting  with  the  Salmon  pulse  at  166 
meters,  is  shown  in  Figure  6.  The  amplitudes  for  the  main  peaks  now  are  in  substan¬ 
tial  agreement  with  the  data  and  the  length  and  amplitude  of  the  precursor  are  also 
reproduced  fairly  well.  Still  there  remains  a  very  abrupt  transition  from  precursor  to 
main  pulse  which  is  clearly  sharper  than  the  experimental  data. 

One  further  refinement  has  been  applied  to  the  partial  shear  failure  model 
in  order  to  avoid  the  abrupt  transition  between  precursor  and  pulse.  Since  there  is 
certain  to  be  a  range  of  material  properties  even  in  the  relatively  homogeneous  salt 
dome,  it  is  reasonable  to  require  a  range  of  thresholds  for  the  shear  failure.  The  model 
has  been  altered  to  produce  a  variation  of  failure  threshold  values  of  compressional 
strain  over  a  range  of  30%  with  a  constant  probability  about  the  10"^  value.  Each  cell 
in  the  finite  difference  calculation  is  given  its  own  threshold  which  is  randomly  drawn 
on  this  basis.  The  set  of  pulses  at  the  Salmon  instrument  ranges  then  calculated  is 
given  in  Figure  7.  The  result  is  a  smoother  transition  from  precursor  to  main  pulse  in 
a  manner  which  is  quite  similar  to  the  actual  Salmon  data  shown  in  Figure  1.  While 
it  is  probably  possible  to  achieve  a  detailed  fit  to  the  data  by  further  such  refinements, 
this  is  not  a  very  meaningful  thing  to  do  since  the  mechanisms  are  not  understood  to 
the  required  level  of  detail.  The  important  point  is  that  it  is  possible  to  reproduce 
the  data  to  a  substantial  degree  using  only  a  few  physically  based  parameters. 
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Figure  6.  Calculated  pulses  at  Salmon  ranges  for  shear  failure  and  Q=10. 
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Figure  7.  Calculated  pulses  at  Salmon  ranges  for  smoothed  shear  failure 
and  Q=10. 


SECTION  4 


COMMENTS 


This  work  indicates  that  the  elastic  precursor  or  leading  toe  seen  in  Salmon 
near-field,  moderate  strain,  velocity  data  is  reproduced  rather  well  with  the  hypothesis 
of  partial  shear  failure  which  is  activated  for  the  duration  of  the  pulse  when  the 
compressional  strain  exceeds  10"^.  This  also  gives  an  attenuation  mechanism  which 
accounts  for  much  of  the  energy  loss  seen  in  the  decay  of  pulses  from  Salmon  with 
range.  However,  the  overall  attenuation  produced  is  not  quite  adequate  to  account  for 
that  seen  in  the  data.  The  addition  of  a  linear  absorption  band  attenuation,  which 
is  active  over  much  of  the  significant  frequency  range  appropriate  to  Salmon  and 
which  has  a  Q  of  10,  then  provides  a  propagation  model  which  nearly  reproduces  the 
signals  at  ranges  beyond  166  meters  when  the  observed  signal  at  this  range  is  used 
as  the  source.  Furthermore,  this  threshold  mechanism  provides  a  transition  to  more 
modest  attenuation  at  small  strains  which  is  required  to  be  consistent  with  Cowboy 
data.  While  there  is  no  assurance  that  this  mechanism  applies  in  other  than  the  salt 
medium.  Ferret  points  out  that  “elastic  precursors”  of  a  similar  character  have  been 
seen  in  UGT  pulses  in  both  alluvium  and  dolomite. 

In  addition  to  accounting  for  much  of  the  Salmon  data,  the  shear  failure 
mechanism  also,  when  applied  to  different  yield  events  in  salt,  will  produce  simple 
scaling  as  observed  over  a  wide  range  of  explosive  events. 

The  fact  that  the  reduction  in  compressional  wave  speed  is  attributed  to 
shear  failure,  rather  than  alteration  of  some  other  material  property,  is  largely  a 
matter  of  consistency  with  past  thinking  on  modes  of  material  behavior;  there  is  no 
direct  experimental  link  to  shear  properties.  The  general  agreement  with  data  which 
results  here  could  just  as  well  have  been  produced  by  any  method  which  reduces  the 
compressional  modulus  in  the  required  amount. 

The  fact  that  the  NER  experiments  give  an  apparently  linear  Q  of  about  20 
for  small  strains  of  10~^  while  Sterling  data  at  comparable  strains  have  much  larger 
Q  is  a  matter  of  some  concern.  However  it  known  that  the  attenuation  of  a  material 
can  be  a  strong  function  of  the  liquid  content.  Spencer’°  made  a  series  of  attenuation 
measurements  in  several  types  of  rock  in  dry  and  water  saturated  states.  He  found  that 
dry  rocks  tend  to  show  only  modest  attenuation  while  the  water  saturated  samples 
had  minimum  Q’s  of  from  4  in  sandstone  to  about  40  in  limestone.  He  found  that  the 
attenuation  for  saturated  rocks  is  quite  frequency  dependent  with  peaks  at  10  to  1000 
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Hz  depending  on  the  material.  He  further  found  that  the  attenuation  spectrum  is 
consistent  with  dispersion  in  the  modulus  so  that  the  mechanism  is  apparently  linear. 
This  sort  of  result  indicates  that  much  care  must  be  taken  in  comparing  experimental 
results  taken  for  different  levels  of  water  saturation  and  at  different  frequencies  even 
when  strain  levels  are  quite  low. 

There  is  a  second  area  where  different  experiments  appear  to  yield  different 
results.  The  work  of  Tittmann  on  decaying  oscillations  of  salt  and  other  bars  give 
apparent  nonlinear  behavior  but  they  also  show  vary  mild  attenuation  in  the  linear 
regime.  As  observed  by  Coyner,  the  attenuation  from  these  multiple  cycle  experiments 
seem  to  be  consistently  much  less  that  from  propagation  (single  pulse)  experiments  at 
comparable  strains  and  frequencies  or  from  hysteresis  loop  experiments.  For  nonlinear 
attenuation,  th«fe  does  not  need  to  be  any  particular  relation  between  the  effective  Q 
from  pulse  and  oscillatory  experiments  -  the  detailed  nature  of  the  mechanism  must 
be  known  to  relate  these.  For  the  example  of  partial  shear  failure,  multiple  cycle 
experiments  will  produce  very  small  attenuation  estimates  since  all  the  failure  will 
take  place  by  the  time  a  single  cycle  has  been  completed  -  after  that  only  the  intrinsic 
attenuation  will  be  active. 

It  should  be  emphasized  that  the  use  of  a  Q  description  for  attenuation 
which  is  clearly  nonlinear  is  probably  not  useful,  except  as  a  crude  estimate  of  the 
degree  of  amplitude  reduction.  In  fact,  the  use  of  Q  in  nonlinear  cases  can  provide  very 
deceptive  results  since  the  frequency  and  range  dependence  are  generally  mixed.  It 
would  be  best  to  describe  the  attenuation  mechanism  directly  in  terms  of  an  equation 
of  state.  This  will  allow  the  application  to  different  experiments.  Put  another  way,  the 
result  of  a  nonlinear  attenuation  for  a  propagation  experiment  can  always  be  used  to 
find  an  effective  Q  by,  say,  use  of  spectral  ratios  but  the  same  mechanism  will  produce 
a  different  effective  Q  if  the  character  of  the  initial  waveform  is  changed.  The  effective 
Q  is  not  a  robust  quantity  for  nonlinear  attenuation. 

We  expect  to  use  the  finite  difference  propagation  code,  which  was  developed 
as  a  tool  for  our  current  work  on  Salmon,  as  a  testbed  for  development  of  nonlinear 
constitutive  relations.  This  avoids  the  use  of  Q  and  concentrates  directly  on  the 
physical  mechanisms  which  produce  attenuation.  Ideally  one  would  like  to  express  a 
constitutive  relation  in  terms  of  parameters  known  to  be  important,  such  as  porosity, 
effective  stress,  crack  density,  crack  growth,  etc.,  in  order  to  allow  a  model  which  can 
be  generalized  to  different  media. 
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APPENDIX  A 


WILKINS  FINITE  DIFFERENCE  EQUATIONS 


The  difference  method  used  to  solve  for  the  spherical  propagation  of  a  pulse 
is  based  on  that  given  by  Wilkins. ‘  The  essential  features  of  this  approach  are  as 
follows: 

The  equation  of  motion  is 

^  Er  -  So 

V  dr  ^  ^  r 

where  the  principal  stresses  are 

Er  =  -(P  +  q)+s, 

=  —  (P  +  q)  +  sj 

in  terms  of  the  hydrostatic  pressure,  P,  the  artificial  viscous  pressure,  q,  and  the  stress 
deviators.  The  reference  density  is  the  relative  volume  is  V  and  U  is  the  particle 
velocity. 


(A  -  1) 

(A-2) 

(A-3) 


The  continuity  equation  is 

V  1  dCr^U) 

V  “  r2  ar 

where  r  is  spherical  radial  position. 

The  energy  equation  is 


(A -4) 


E  -  V[sifi  +  2S2C2I  +  (P  +  q)V  =  0 
where  E  is  the  internal  energy  and  f  are  strains: 

€1  =  dU/dr 
£2  =  U/r 

Artificial  viscous  pressure  is  taken  from  ViscelH’s’  work  as: 


(A -5) 

(A-6) 

(A-7) 


A-1 


„  A  au 


(A -8) 


where  Cl  is  a  constant  (taken  as  0.5),  a  is  the  sound  speed  and  Ar  is  the  cell  spacing. 


For  the  elastic  case  the  equation  of  state  is 


=  S^-Ui-jV/V) 

i,  =  2M(€,-iv/V) 
P  =  -(^  +  5A‘)V/V 


(A-9) 

(A-10) 

(A-11) 


where  A  and  /i  are  the  Lame  elastic  constants  (A  +  |/x  is  the  bulk  modulus). 


The  finite  difference  version  equations  is  taken  by  division  of  the  material 
into  N  mass  intervals 


3V0 


j  =  1,...,N 


(A  -  12) 


The  equation  of  motion  for  the  velocity  for  the  element  at  the  n***  time 


where 


I  At**  r  1 

UJ>«  =  U,"->  +  —  -  (E,),".  j]  +  2At"(^;)  (A  -  13) 

+  (A.14) 

(EDhs  =  +  (A-**) 

"  i  h  (^) 


(A-14) 

(A-15) 


(A-16) 


A-2 


At  an  outside  regional  boundary  J 


(A  -  26) 


Xj+i  = 


(Atj+5) 

12 


[(<.*)” -(ur*)' 


The  strains  are  defined  by 


t"+1 


"+a  _  "+* 

urV + ur* 

r",V+rr- 


('-Cs 


= 


(A-27) 

(A-28) 


and  the  stress  deviators  are  advanced  by 


n+1  _ 


n+1  _ 


(si)j"+i  +  2/x 

(s2)JVl  +  2/* 


'■+!  3  \  V"+i 

^  *^j+i  3  V  V"+5  /  1 


(A-29) 

(A-30) 


Artificial  viscous  pressure  is  included  by 


n+  Q  _  Q  n+  A 

qj+Z  =  Clop  r;.^/ 


c*  -  ur* 


(A  -  31) 


where  Cl  «  0.5  (constant)  and  a  in  the  compressional  wave  speed. 
The  energy  equation  is  advanced  by 


(E,)r,j +v-^(s, )":*(.■, );;Mt 


n+i 


(E»)'.i 


n+1  _ 


n+1  _ 


(A-32) 

At"*!  (a-33) 

{e;  -  {j|P"*‘  +P”I  +  q"*!}  ■  |V"*'  -  V"|  +  AE,  +  AEj)  ^A-34) 


A-4 


For  the  current  application,  this  energy  is  just  a  diagnostic  and  is  not  fed 
back  through  any  equation  of  state. 

The  hydrostatic  pressure  for  the  elastic  case  is 

5")  (’’"}-») 

REFERENCES 

1.  Wilkins,  M.L.,  “Calculation  of  Elastic-Plastic  Flow,”  Methods  in  Computational 
Physics,  Vol.  3,  1964. 

2.  Viecelli,  J.A.,  “The  Linear  Q  and  the  Calculation  of  Decaying  Spherical  Shocks 
in  Solids,”  J.  Comp.  Phys.,  12,  187  (1973). 


A-5 


APPENDIX  B 


INCLUSION  OF  DAY-MINSTER  Q  IN  FINITE  DIFFERENCES 


Day  and  Minster’  have  shown  how  to  include  an  arbitrary  linear  Q  function 
in  time  stepping  calculations.  For  an  absorption  band  an  analytic  solution  is  available. 
.\n  outline  of  their  methods,  as  easily  generalized  to  the  spherical  case  is  given  here. 

For  a  single  normalized  relaxation  function,  m(t),  the  stresses  and  strains 
are  related  by 


S.  = 

(A  +  2p)  J  Tn{t  —  T)ti(r)dT  +  2A  J  Tn{i  ~  T)f2(7-)dr 

(B-l) 

2,  = 

(2A  4-  2p)  j  m(t  -  r)c2(’')dr  +  ^  j  *Ti(t  —  t)ti(r)dr 

(B-2) 

Here  the  strains  are 

du 

(B-3) 

-  Tt 

u 

C2  =  - 

r 

(B-l) 

where  u  is  displacement.  Generally  one  could  have  two  different  relaxation  functions 
(e.g.,  bulk  and  shear)  but  we  shall  not  consider  this  possibility.  The  A  and  n  are  the 
usual  Lame  constants  which  are  now  the  unrelaxed  or  high  frequency  moduli  of  the 
medium.  It  stresses  can  be  written  in  terms  of  Q  corrected  strains,  e,  as 


Lr  —  (A  -f-  2p)ei  +  2Ae2  (B-5) 

2}  =  (2A  +  2p)e2  +  Acj  {B*6) 

Day  and  Minster  have  shown  how  to  express  the  ej  in  terms  of  the  cj 

Ci  =  J  m(t  -  r)eidT  (B  -  7) 

using  a  sequence  of  m  Fade  approximants  to  write  the  integral  relation  as  a  differential 
relation.  For  an  absorption  band  attenuation  with  relaxation  times  between  Ti  and  T2 
and  with  a  flat  spectrum  they  show  that  the  integral  relation  can  be  replaced  by 


B-1 


(B  -  8) 


m 

ei(t)  =  «i(t)  -  ^ 


k=l 


where  ^ 

^6*'  +  -~--W„Qo^^  £i(t)  k-1, ...,  m  (B  -  9) 

Here  Qo  is  the  target  Q  in  the  absorption  band,  the 

^  [4(rf '  -  T^')  +  (rf »  +  r,-^)]  (B  -  10) 

where  the  t\i  are  the  abscissas  amd  Wk  are  the  weights  for  m-point  Gauss-Legendre 
quadrature.  The  index  m  is  that  of  the  Fade  approximant.  As  m  increases,  the 
solution  converges  to  the  analytic  result  which  in  the  frequency  domain  is 


•m  ^ 

iu;in(a;)  =  m(a/)  =  1 - 

fQo 


( 


-  n) 

1  +  w*rir2 


The  convergence  of  the  sequence  with  increasing  m  is  shown  in  Figure  B-1.  For 
calculations  given  in  the  main  text  m=5  was  used  since  convergence  is  near.  Use  of 
significantly  larger  m  slows  down  the  calculations  of  pulse  propagation  since  each  m 
adds  an  additional  differential  equation  which  must  be  carried.  In  finite  difference 
form  the  g  are  advanced  for  both  i  by 


At" 


Tj  -  r, 


2  +  At"t\  xQo 


W 


‘)  (v1 


+  V})  (B-12) 


and  then  the 


e  = 


k=l 


(B  - 13) 


These  are  then  used  to  advance  the  stresses  by 


B-2 


so 


Si  =  ~  ca)  (B-14) 

P  =  (^  +  |»‘)(ei  +  %)  (B-15) 

=  («.)?*  (B.16) 

P|”4  =  P^}  +  (B-") 

(..),■;!  =  (b-18) 


are  used  in  place  of  the  procedure  applied  in  Appendix  A. 
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